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SUMMARY 

Analysis  of  the  performance  characteristics  of  a labo- 
ratory model  real-time  thermoplastic  optical  data  modulator 
was  accomplished.  Design  of  the  modulator  was  performed  in 
a previous  program.  Experimental  data  for  bandwidth,  dynamic 
range,  spatial  resolution  and  optical  quality  were  obtained 
indicating  a 1evel  of  performance  suitable  for  entry  of  wide- 
band, high  dynamic  range  data  into  a coherent  optical  pro- 
cessor or  optical  computer,  or  use  of  the  modulator  as  a large 
screen  display  device. 

The  thermoplastic  modulator  provides  electron  beam  writing 
on  a reusable  deformable  thermoplastic  surface  giving  a two- 
dimensional  optical  phase  recording  which  may  be  rapidly  erased 
upon  command.  Data  written  on  the  plastic  surface  may  be  re- 
tained for  arbitrary  lengths  of  time  down  to  a minimum  of  about 
100  milliseconds,  prior  to  erasure.  A 70  MHz  temporal  band- 
width, 30  c/mm  spatial  bandwidth  and  a recording  dynamic  range 
of  at  least  40  db  were  realized.  The  experimental  modulator 
system  was  designed  in  two  configurations  which  differed  mainly 
in  the  thermoplastic  recording  surface  geometry.  These  two 
implementations  were  a fixed  plate  and  a rotated  disc  recording 
arrangement.  The  fixed  plate  system  provides  a 5 inch  diameter 
substrate  on  which  central  3 inch  diameter  area  is  used  for  re- 
cording. The  rotated  disc  system  utilizes  a 15  inch  diameter 
substrate  which  is  mechanically  rotated  past  the  repetitively 
line  scanned  electron  beam.  The  thermoplastic  recording  region 
lies  in  a band  between  radii  of  2 inches  and  6.5  inches.  Various 
data  formats  are  provided  by  electronic  scan  control  of  the 
writing  electron  beam  in  addition  to  the  mechanical  rotation 
of  the  recording  surface  in  the  case  of  the  rotated  disc  system. 
Optical  quality  of  the  recording  medium  of  better  than  one-fourth 
wavelength  flatness  over  a one  inch  region  is  achieved  using 
materials  and  thermoplastic  coating  techniques  developed  in  this 
program. 
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A thermoplastic  modulator  which  provides  real-time  inser- 
tion of  data  to  the  input  of  a coherent  optical  processing 
channel  has  been  designed  and  its  performance  properties  eva- 
luated. The  initial  fabrication  of  components  for  this  experi- 
mental system  was  accomplished  in  a previous  program^ 

Coherent  optical  data  processing  systems  offer  excellent  and 
often  unique  performance  features.  The  use  of  coherent  optical 
systems  of  course  requires  entry  of  the  data  to  be  processed  in 
a two-dimensional  spatial  format  at  the  input  plane  of  the  pro- 
cessing optics.  Where  the  signal  to  be  processed  is  a time  vary- 
ing voltage,  a recording  device  (data  modulator  or  transducer)  is 
required  that  will,  generally,  be  comprised  of  a writing  beam  that 
is  modulated  in  proportion  to  the  signal  to  be  recorded,  and,  a 
recording  medium  on  which  the  beam  writes  in  the  desired  spatial 
format.  For  optical  processing  the  recorded  data  is  illuminated 
with  a coherent  light  beam,  the  recording  causing  a variation  in 
either  the  phase  or  amplitude  of  the  illuminating  beam  depending 
upon  the  type  of  recording  medium  employed.  This  spatially  modu- 
lated light  beam  is  then  processed  using  lenses  and  spatial  fil- 
ters to  provide  the  data  processing  or  display  functions  of  in- 
terest. 

There  are  several  concepts  potentially  available  for  appli- 
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cation  as  a near  real-time  data  recorder  having  high  optical 
quality  and  a reusable  recording  medium.  The  recording  con- 
cepts may  be  categorized  generally  as  to  type  of  writing  beam 
and  recording  material  used.  Typically  either  an  electron  beam 
or  light  beam  is  used  for  writing.  The  list  of  recording  medium 
. candidates,  which  is  somewhat  dependent  on  the  type  of  writing 

beam  used,  includes  deformable  plastics,  deformable  membranes, 
variable  scattering  in  liquid  crystals  or  ferroelectric  ceramics, 
variation  of  optical  transmittance  in  cathodochromics  or  photo- 


chromics,  and,  variation  of  induced  birefringence  in  several 
possible  electro-optic  materials. 

The  use  of  a thermoplastic  data  modulator  employing  an  elec- 
tron writing  beam  offers  considerable  flexibility  in  design  of 
the  beam  scanning  provisions  particularly  for  the  case  of  a 
two-dimensional  scan  format.  In  addition,  modulation  of  the 
writing  beam  in  correspondence  with  the  input  signal  to  be 
recorded  is  comparatively  straightforward  for  bandwidths  to  at 
least  one  hundred  megahertz  and  can  be  accomplished  to  consider- 
ably greater  bandwidths  with  special  electron  beam  electrode 
design.  These  factors  provide  a natural  motivation  for  use  of 
a real-time  modulator  which  employs  an  electron  beam.  The 
suitability  of  the  thermoplastic  recording  material  is  governed 
by  such  performance  factors  as  recording  quality  and  sensitivity, 
optical  quality,  optical  diffraction  efficiency  and  recycle  life. 
Thermoplastic  offers  the  potential  for  good  recording  sensitivity 
and  quality  with  optical  readout  allowable  over  a wide  range  of 
wavelengths  and  at  comparatively  good  efficiency.  A deformable 
thermoplastic  has  been  developed  as  a near  real-time  recording 
medium  with  good  optical  quality  and  thickness  uniformity.  Data 
recording  is  accomplished  by  causing  surface  deformations  of  a 
thin  thermoplastic  layer  as  depicted  in  Figure  1.  The  deforma- 
tions are  caused  by  the  electrostatic  force  of  the  charge  de- 
posited by  the  writing  electron  beam  when  the  plastic  is  suf- 
ficiently softened  by  raising  its  temperature.  \ further  in- 
crease in  temperature  causes  the  recorded  data  to  erase  by 
allowing  the  surface  charge  to  be  conducted  through  the  plastic 
with  surface  tension  forces  then  acting  to  smooth  the  surface. 

The  readout  laser  beam  becomes  phase  modulated  in  correspondence 
with  the  recorded  variations  in  thickness  of  the  deformed  plas- 
tic as  depicted  in  the  sketch  of  Figure  2. 
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FIGURE  1.  DATA  STORAGE  BY  THERMOPLASTIC  DEFORMATION 


FIGURE  2. 


DATA  READOUT  FOR  OPTICAL  PROCESSING 
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SYSTEM  CONFIGURATION 


The  thermoplastic  data  modulator  developed  is  designed 
in  two  configurations  which  differ  mainly  in  the  recording 
plate  geometry  as  shown  in  Figures  3 and  4.  Figure  3 shows 
a configuration  in  which  the  TP  recording  plate  is  fixed  in 

j 

position  much  as  the  phosphor  faceplate  of  a cathode  ray  tube. 

A larger  TP  recording  plate  which  can  be  mechanically  rotated 
either  continuously  or  incrementally  is  used  in  the  second  con- 
figuration shown  in  Figure  4.  At  present  the  fixed  plate  system 
makes  use  of  a 5 inch  diameter  fused  silica  substrate  with  a 
transparent  gold  coating  over  which  the  thermoplastic  is  coated 
in  the  central  4 inch  diameter  region.  In  the  rotated  plate  sys- 
tem a 15  inch  diameter  substrate  with  a thermoplastic  layer  in  a 
radial  interval  between  radii  of  2 and  6.5  inches  is  used. 

The  system  consists  generally  of  an  evacuated  envelope  (con- 
taining the  electron  gun,  the  TP  recording  plate,  and  a mirror) 
and  supporting  electronics  and  vacuum  pumping  provisions.  The 
expanded  laser  beam  used  for  data  readout  is  passed  through  an 
entrance  window  of  the  evacuated  recording  chamber  and  reflected 
from  the  internal  mirror  so  as  to  pass  through  the  TP  recording 
plate. 

The  electron  beam  is  generated  with  a triode  electrode  struc- 
ture in  which  a replaceable  indirectly  heated  dispenser  cathode  is 
used  and  beam  focus  is  accomplished  with  electrostatic  lenses. 

Beam  current  is  modulated  by  varying  the  grid  to  cathode  voltage. 

Beam  scanning  is  obtained  with  an  orthogonal  set  of  magnetic  deflec- 
tion coils.  A 6 KV  beam  voltage  with  beam  current  up  to  about 
lya  is  used. 

In  the  fixed  plate  system  the  writing  beam  can  be  electroni- 
cally scanned  in  a rectangular  raster  or  a polar  raster  format. 

With  the  reotating  disc  configuration  the  electron  beam  scans 
repetitively  over  a single  line  as  the  disc  rotates  past  it  or 
in  a sequence  of  small  rasters. 
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FIGURE  3.  THERMOPLASTIC  MODULATOR  - FIXED  PLATE 
CONFIGURATION 


FIGURE  4.  THERMOPLASTIC  MODULATOR  - ROTATING  DISC 
CONFIGURATION 
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The  storage  medium  is  comprised  of  a thin  (-6  ym)  thermo- 
plastic layer  coated  onto  a glass  substrate,  with  a thin  trans- 
parent conductive  coating  between  as  shown  in  Figure  1.  The 
conductive  film  serves  as  a plane  of  constant  electrical  poten 
tial.  It  also  serves  as  a means  for  electrical  heating  (I  R)  of 
the  thermoplastic  for  rapid  deformation  development  and  data 
erasure.  A voltage  is  applied  across  opposite  edges  of  the  con- 
ductive coating  for  heating  purposes. 

Thermoplastic  can  be  operated  either  in  a continuous  develop- 
erase  mode  or  a pulsed  develop -erase  mode.  In  the  continuous 
mode  the  TP  temperature  is  held  at  a level  which  allows  fairly 
rapid  development  of  recorded  deformations  and  slightly  slower 
erasure  of  the  recorded  data.  Some  variation  of  develop  and 
erase  time  constants  is  possible  depending  upon  the  material 
used  and  the  operating  temperature  selected.  For  a temperature 
of  60°C  the  system  described  here  provides  deformation  development 
in  about  5 seconds  and  a 60  second  erasure.  Pulsed  control  of 
the  TP  temperature  allows  development  in  less  than  one  second 
after  which  the  recorded  data  can  be  retained  until  application 
of  an  erase  pulse  which  will  rapidly  remove  the  stored  data. 

Development  and  erasure  of  signals  written  on  the  thermo- 
plastic are  governed  by  the  physical  and  electrical  properties 
of  the  TP  material  and  their  change  with  temperature.  The  charge 
deposited  on  the  thermoplastic  surface  during  writing  results  m 
an  electric  field  force  between  the  TP  surface  and  the  equipoten- 
tial  plane  at  the  conductive  coating.  This  force  acts  to  deform 
the  surface.  Raising  the  thermoplastic  temperature  allows  it  to 
change  from  a solid  to  a rubber-like  condition  and  then,  with  fur- 
ther temperature  increase,  to  a liquid  state.  Deformation  devel- 
opment can  occur  while  the  thermoplastic  is  in  either  of  these 
softened  states.  The  depth  of  the  deformation  will  be  limited  by 
counteractive  surface  tension  and  viscoelastic  forces.  Erasure 
of  signal  deformations  occurs  when  the  surface  charge  is  removed 
and  the  thermoplastic  temperature  is  elevated  so  the  plastic 
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reaches  a non-solid  state  and  can  be  re-distributed  so  as  to 
fc  o a smooth  surface  in  response  to  the  surface  tension  and 
viscoelastic  restoring  forces.  In  erasure,  the  written  surface 
charge  is  removed  by  electrical  conduction  through  the  plastic 
to  the  conductive  layer  when  the  TP  resistance  is  lowered  by 
heating. 


DATA  READOUT 

Data  written  on  the  thermoplastic  is  read  out  with  a laser 
beam  expanded  to  illuminate  the  area  over  which  stored  signals 
of  interest  appear,  as  shown  in  Figure  2.  The  sketch  shows 
readout  by  transmission.  (Reflective  methods  are  also  possible). 
The  basic  properties  of  optical  readout  of  the  phase  data  is  dis- 
cussed in  this  section  with  a more  detailed  description  relating 
to  recording  format  given  later  in  Section  6. 

In  rectangular  coordinates  x,  y the  complex  transmittance  of 
the  thermoplastic  plate  can  be  written  as  T (x  ,y)  e^  a where  T 

is  the  amplitude  transmittance  and  a the  phase  shift  caused  by 
the  plastic  and  substrate.  Given  an  illuminating  beam  having  a 
wavefront  with  uniform  phase  and  constant  amplitude  over  an  aper- 
ture A(x,y),  the  emergent  beam  S(x,y)  at  the  output  side  of  the 
glass  substrate  will  be  of  the  form 


S(x,y)  = A(x,y)  T(x, y)e^ a ^x' 


The  amplitude  transmittance  will  normally  be  constant,  i.e., 

T(x,y)  = Tq.  Phase  shift  at  the  optical  wavelength  can  be 
written  as  a(x,y)  = [2TT(n  - 1)/A]d(x,y)  + 3,  where  d is  surface 
deformation  depth,  n is  the  index  of  refraction  of  the  thermoplastic 
and  3 is  a phase  shift  caused  by  materials  other  than  the  thermo- 
plastic, i.e.,  substrate  and  conductive  coating.  Thus  the  emer- 


For  example,  consider  the  deformation  d to  be  a one -dimensional 
sinusoid  of  spatial  frequency  w 


d = D + Dsin  w x 
o x 


as  indicated  in  Figure  2.  Considering  the  x variation  only  for 
convenience,  the  expression  for  the  phase  modulated  readout  beam 
S can  be  written  as 
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where  we  have  ignored  constant  phase  terms  and  J^[  ] is  a Bessel 
function  of  the  first  kind,  order  m.  An  example  of  the  amplitude 
variation  given  by  the  Bessel  functions  is  given  in  the  data  of 
Figure  5.  The  indexing  variable  m takes  on  all  integer  values. 
Thus,  there  are  plane  waves  e^maixX  emanating  from  the  modulator; 
each  has  a field  amplitude  J 


m 


tir, 

2 • (n  - Jd|  and  all  come  through 


x r\ 

the  aperture  A(x) . 

If  the  transducer  output  is  Fourier  transformed  optically  as 
depicted  in  Figure  2,  the  Fourier  transform  lens  will  focus  each 
of  the  plane  waves  e^maixX  into  a small  spot  of  light  at  the  output 
transform  plane  of  the  lens . In  this  example  each  spot  of  light 
will  have  a shape  or  spatial  distribution  dependent  on  the  size 
and  shape  of  the  data  aperture  or  apodization  within  the  aper- 
ture A(x) . For  a uniformlv  i lluminated  round  aperture  the  familiar 


Airy  disc  distribution, 


Jlt  1 


i ] 


is  obtained.  The  peak  value  of 


the  light  int  msity  of  the  focused  spot  will  be  proportional  to 
J^{[2-rr(n  - 1) / X ]D)  and  thus  depends  upon  the  depth  D of  the  sinu- 
soidal deformation  in  the  thermoplastic.  The  distance  of  a spot 
from  the  zero  order  position  will  be  directly  proportional  to  the 
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FIGURE  5.  BESSEL  FUNCTIONS  OF  THE  FIRST  KIND 
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spatial  frequency  in  the  plane  wave  expression  e-'nuxX  and  its 

angular  location  will  be  on  a line  orthogonal  to  the  input  grating 

lines.  Only  data  corresponding  to  one  of  the  first  orders  (m  = +1) 

are  usually  used  for  optical  processing.  The  deformation  depth  D 

is  then  chosen  to  optimize  the  performance  of  the  transducer  in 

terms  of  linearity  and  diffraction  efficiency.  Diffraction  effi- 

2 

ciency  for  first  order  data  is  taken  as  the  ratio  J,  [ ] to 

2 J f ] . 
m m1  J 

While  phase  modulation  of  the  readout  beam  provides  greater 
diffraction  efficiency  than  other  modulator  concepts  it  will 
typically  be  held  to  a few  per  cent  to  obtain  linear  operation 
and  to  limit  higher  order  diffraction  and  intermodulation  products. 
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Several  system  as  well  as  component  characteristics  have 
been  investigated  and  some  of  the  results  will  be  reviewed  in 
this  section. 

System  performance  was  tested  while  using  the  fixed  plate 
modulator  in  an  optical  Fourier  transform  data  processing  mode 
of  operation  as  depicted  in  Figure  6.  Data  recorded  as  thermo- 
plastic surface  deformations  was  illuminated  with  an  expanded 
and  collimated  HeNe  laser  beam.  The  phase  modulated  light  beam 
wavefront  emergent  from  the  TP  layer  is  Fourier  transformed  by 
lenses  as  shown  in  Figure  6.  The  output  light  distribution  at 
the  spatial  frequency  (Fourier  transform)  plane  was  examined 
as  a function  of  the  input  voltage  used  to  modulate  the  current 
of  the  electron  beam  and  several  operating  parameters.  Beam 
current  was  modulated  relative  to  a bias  level  set  by  a bias 
voltage  at  the  modulation  electrodes  of  the  electron  gun. 
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The  amplitude  modulation  transfer  function  (AMTF)  is  described  by 
the  data  of  Figure  7.  Shown  is  the  relative  amplitude  of  the  peak 
value  of  the  first  order  diffracted  spot  of  light  in  the  spatial  fre- 
quency plane  as  a function  of  the  frequency  of  the  input  signal  used 
to  modulate  writing  beam  current.  Light  amplitude  was  computed  as  the 
square  root  of  measured  light  intensity.  The  data  was  taken  at  approxi- 
mately 2,  5 and  10  MHz  increments  for  curves  A,  B and  C,  respectively. 
The  AMTF  curves  corresponding  to  several  different  values  of  writing 
beam  sweep  speed  are  given.  Though  not  shown  in  the  data,  a recorded 
response  was  realizable  up  to  frequencies  of  90  MHz  when  using  the 
highest  sweep  speed.  A temporal  bandwidth  of  about  70  MHz  at  the  5 
per  cent  response  points  is  achieved  when  using  the  highest  speed.  As 
sweep  speed  is  reduced  the  AMTF  is  scaled  toward  lower  frequencies 
and  at  a higher  sweep  speed  a greater  bandwidth  can  be  realized  to 
the  extent  that  sufficient  beam  current  is  available  for  exposure  of 
the  thermoplastic. 

Relative  gain  and  dynamic  range  data  for  the  system  can  be  seen 
in  Figure  8.  System  gain  is  given  by  the  plot  of  relative  amplitude 
of  the  peak  value  of  the  first  order  diffracted  light  spot  in  the 
spatial  frequency  plane  as  a funtion  of  the  peak-to-peak  value  of 
the  input  sinusoidal  voltage  used  to  modulate  the  writing  electron 
beam.  Frequency  of  the  input  voltage  was  held  constant  at  either 
10  MHz  or  20  MHz  when  performing  this  test.  The  amplitude  of  the 
peak  value  for  the  noise  background  and  of  the  second  order  dif- 
fracted light  at  the  spatial  frequency  plane  is  also  given.  As 
shown  a recorded  signal  variation  for  inputs  from  0.1  to  at  least 
18  volts  is  achievable.  The  ratio  of  maximum  oucput  (22)  to  RMS 
noise  (0.2)  for  the  modulator  was  110  or  a dynamic  range  of 
40 . 8 db  . 

Diffraction  efficiency,  taken  as  the  ratio  of  first  order  to 
zero  order  diffracted  light  intensity,  is  about  3 per  cent  for  the 
recording  of  a 15  volt  peak-to-peak  input  sinusoidal  voltage. 

An  indication  of  the  quality  of  the  diffracted  light  spot  in 
the  spatial  frequency  plane  is  shown  by  the  photographs  of  Figure  9. 

The  light  distribution  shown  in  Figure  9a  and  9b  is  caused  by  a 30 
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MHz  recording  (sine  wave  grating)  within  a square  illumination 
aperture  at  the  thermoplastic  surface.  Figures  9c  and  9d  are  dif- 
fraction patterns  obtained  with  no  recorded  data  (zero  frequency) 
with  a square  input  aperture  and  Figures  9e  and  9f  are  similar  data 
obtained  when  using  a circular  input  aperture.  Good  recording 
medium  and  recorded  data  quality  are  indicated  by  the  light  distri- 
butions shown  in  Figure  9 which  correspond  to  nearly  diffraction 
limited  performance  for  data  apertures  up  to  at  least  an  inch. 

A microscopic  scan  of  the  diffraction  pattern  intensity  distri- 
butions for  the  circular  and  square  apertures  is  shown  in  Figures 
10  and  11  which  indicate  good  correspondence  to  theoretically  pre- 
dictable results  in  mainlobe  and  sidelobe  relationships.  The  left 
and  right  sidelobes  are  nonsymmetric  in  peak  levels  which  is  be- 
lieved caused  by  a slight  nonuniformity  in  coating  thickness.  The 
magnified  frequency  plane  spots  were  scanned  with  a photometer  having 
a 10  pm  circular  scanning  aperture.  The  photometer  output  was  ampli- 
fied with  a logrithmic  amplifier  and  recorded  on  an  X-Y  chart  re- 
corder to  obtain  the  data  shown  in  these  figures.  The  vertical 
(intensity)  scale  is  logrithmic  with  approximately  0.37  dB  per 
smallest  division.  The  horizontal  (distance)  scale  is  1.2  pm 
per  smallest  division  for  the  1"  aperture  data  and  about  1.6pm 
per  division  for  the  1/2"  aperture  data  when  observed  at  the 
output  (frequency)  plane  of  a diffraction  limited  12"  focal 
length  lense.  The  intensity  level  baseline  is  caused  by  the 
photometer  noise  level.  The  single  frequency  spot  distributions 
shown  in  these  figures  correspond  to  the  theoretical  diffrac- 
tion limited  distributions  to  within  about  15  percent  indicating 
good  optical  quality  and  writing  fidelity. 
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FIGURE  9.  SINGLE  FREQUENCY  SPOT  DISTRIBUTION  (C-Circular  Aperture 
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FIGURE  10.  MICROSCOPIC  SCAN  OF  FREQUENCY  PLANE  FOR 
CIRCULAR  INPUT  APERTURE 
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FIGURE  II . MICROSCOPIC  SCAN  OF  FREQUENCY  PLANE  FOR  SQUARE 
INPUT  APERTURE 
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THERMOPLASTIC  AND  SUBSTRATE  PROPERTIES 


Development  of  fabrication  methods  for  the  thermoplastic  re- 
cording medium  was  continued  during  this  program  using  only  one 
of  several  potentially  useful  plastic  materials.  The  main  em- 
phasis has  been  on  design  of  the  substrate  and  methods  for  appli- 
cation of  the  plastic  layer  in  the  range  of  5 to  7 ym  thick  with 
a high  uniformity  of  coating  thickness.  Two  substrate  sizes 
have  been  developed:  a 5"  diameter  substrate  on  which  a 4" 

diameter  TP  coating  is  applied  for  the  fixed  plate  system;  and 
a 15"  diameter  substrate  that  is  coated  in  a band  between  radii  of 
2"  and  6.5"  for  the  rotating  disc  system.  These  substrates  are 
constructed  from  3/8  thick  optical  quality  fused  silica  with  a 
transparent  gold  conductive  film  and  silver  electrode  strips  for 
electrical  contact  as  shown  in  Figure  12.  The  coating  methods  de- 
veloped for  the  5"  diameter  substrate  have  been  described  previously^ 
and  only  the  large  disc  is  discussed  in  detail  in  this  section. 

With  the  large  disc  many  of  the  problems  that  occurred  with  the 
previously  developed  5"  diameter  plate  are  magnified  and  other 
factors  that  were  not  a problem  before  became  a factor.  For  a 
satisfactory  coating, good  flatness  and  low  particle  contamination 
are  required.  Let  us  first  consider  producing  the  best  flatness. 

To  improve  flatness  the  substrate  has  to  be  flat  itself.  Due  to 
manufacturing  factors  a substrate  as  flat  as  the  smaller  discs  is 
difficult  to  obtain.  A precision  granite  flat,  Figure  13„  was  im- 
plemented as  the  substrate  support.  This  not  only  offered  overall 
support  for  the  substrate,  but  also  provided  a flat  surface  for 
the  glass  to  conform  to  in  cases  in  which  the  substrate  was  slightly 
warped  in  some  fashion.  The  substrate  can  be  leveled  to  at  least  10 
seconds  of  arc. 

As  with  the  smaller  coatings  a certain  solvent  and  solvent  layer 
thickness  are  required  with  a given  thermoplastic.  These  parameters 
may  vary  with  thermoplastic,  but  appear  to  be  independent  of  coating 
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FIGURE  12.  THERMOPLASTIC  SUBSTRATES 


FIGURE  13.  COATING  TABLE  AND  FORMS  (The  outer  form  is  removed  showing 
the  inner  retaining  forms  that  define  the  coating  layer 
boundary  on  the  substrate.  In  the  background  is  the  front 
of  the  clean  room  from  which  clean  laminar  air  flows  hori- 
zontally. To  the  left,  not  shown,  is  the  vapor  hood, 
which  collects  the  escaping  solvent  vapors.  The  form 
is  sitting  on  a precision  granite  plate.  The  plate  is 
leveled  by  four  screws  (40  threads/inch)  at  the  bottom 
of  the  support  stands). 


area.  The  thermoplastic  solution  used,  due  to  a property  of  the 
thermoplastic,  did  not  read  .ly  wet  the  surface  but  had  to  be 
spread  over  the  area  with  a stir  rod,  Figure  14.  First,  the  solu- 
tion was  applied  more  or  less  uniformly  about  the  circumference 
of  the  annulus  being  coated.  A stir  rod  was  used  to  move  the 
solution  layer  to  the  boundaries.  This  joining  was  done  last  in 
order  to  minimize  the  quantity  of  solution  captured  by  the  edge. 

The  spreading  of  the  solution  by  the  stir  rod  required  a technique 
that  applies  the  solution  rapidly  but  uniformly.  It  was  found 
that  the  size  of  the  stir  rod  was  important.  A rod  too  small 
(in  terms  of  the  width  of  the  rod  which  interfaces  with  the  so- 
lution) increased  the  application  time  unduly.  A rod  too  wide 
left  an  area  behind  the  center  of  the  rod  where  the  solution  layer 
was  too  thin.  The  rod  would  appear  to  drag  the  liquid  rather  than 
push  it.  The  rod  should  be  moved  rapidly  around  the  coating  area 
without  overconcern  as  to  how  thick  the  solution  layer  is.  Where 
the  layer  is  too  thin  such  as  occurred  with  the  extra  wide  stir 
rod,  the  solution  layer  will  break  through  after  several  seconds 
so  that  the  area  must  be  wetted  again  from  neighboring  areas. 

Holding  the  solution  layer  in  a quiet  state  after  spreading  and 
before  drying  helps  insure  good  flatness. 

During  the  spreading  period  there  is  a certain  amount  of  sol- 
vent evaporation.  This  is  undesirable  in  that  when  the  solution 
layer  is  in  the  form  of  a large  sessile  drop  the  greatest  evapora- 
tion occurs  near  the  edge,  since  the  vapors  readily  diffuse  to  the 
surrounding  uncoated  area.  Hence,  nonuniformities  in  concentra- 
tion will  be  produced.  Then,  even  with  an  initially  smooth  solution 
layer,  the  final  surface  will  not  be  flat.  To  reduce  the  evapora- 
tion, the  surrounding  atmosphere  is  saturated  with  the  solvent 
vapor.  In  addition,  the  circulation  of  air  is  reduced  by  a double 
enclosure  surrounding  the  solution  layer.  First  is  the  metal 
form  which  also  acts  to  dam  the  solution  layer,  Figure  13,  and  then 
a surrounding  f iberglas-epoxy  enclosure  that  reduces  the  circulation 
into  the  first  form,  Figure  15 . 
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FIGURE  14.  THERMOPLASTIC  APPLICATION  (OUTER  RETAINING  FORM  REMOVED) 
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The  flatness  of  the  solution  layer  can  be  checked  by  ob- 
serving the  reflections  of  the  overhead  lights.  When  it  appears 
that  the  coating  is  sufficiently  flat  so  that  no  more  breakthrough 
will  occur,  the  lids  of  both  forms  are  put  on  for  about  20  minutes. 

Next,  the  coating  is  allowed  about  two  hours  to  dry.  During  that 
period  the  rate  of  circulation  is  gradually  increased.  If  the 
drying  period  is  too  long  a radial  wedge  can  be  introduced  into  the 
coating.  This  occurs  because,  typically,  the  outer  edge  tends 
to  evaporate  more  quickly  so  the  solution  drains  from  inside 
causing  the  solution  layer  on  the  inside  to  become  thin  with  the 
possiblity  of  breakthrough.  The  latter  effect  would  occur  sub- 
sequent to  the  occurrence  of  a wedge  if  it  occurs  at  all,  and  would 
require  a new  start  in  the  procedure.  This  possibility  is  diminished 
by  loading  the  outer  form  with  solvent  in  a trough  along  the  wall 
which  saturated  the  air  with  vapor.  Having  a configuration  in  which 
there  is  even  evaporation  is  essential  but  it  is  hard  to  maintain 
for  a long  time  interval.  The  slower  the  evaporation  and  circulation 
the  more  likely  will  be  a vapor  gradient.  These  problems  were  reduced 
by  the  above  mentioned  addition  of  solvent  vapor  and  were  further 
diminished  by  increasing  the  circulation  gradually  as  the  solution 
layer  becomes  concentrated  and  more  viscous.  The  sequence  used  to 
realize  the  increased  circulation  is  first,  the  lid  of  the  outer 
form  is  placed  ajar;  second,  the  same  lid  removed.  The  latter  must 
be  done  at  the  right  time  and  with  care  in  order  that  no  solvent  is 
dropped  on  the  coating.  The  TP  solution  along  the  edge  should  not 
be  so  dry  that  strings  of  TP  are  pulled  up  and  fall  onto  the  coating, 
or  worse,  that  the  form  is  frozen  onto  the  substrate.  On  the  other 
hand,  it  should  not  be  so  wet  that  the  TP  solution  drips  onto  the 
coating . 

Problems  with  particle  contamination  have  been  considered  and 
much  progress  has  been  made.  The  acetone -acetate  method  whereby 
acetone  is  applied  and  then  an  acetate  sheet  (the  acetate  is  par- 
tially dissolved  and  surrounds  and  traps  the  particles)  which  worked 
on  the  smaller  discs  was  not  as  satisfactory  on  the  large  disc.  With  the 
large  disc,  the  application  time  increased  and  during  that  time  the  ace- 
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i tate  swelled  sufficiently  to  make  a smooth,  bubble  free  acetate 

layer  nearly  impossible.  Changing  to  methyl-ethyl-ketone  improved 
the  results;  however,  the  acetate  would  separate  from  the  substrate 
before  completely  drying.  The  hole  in  the  middle  of  the  large 
disc  substrate  and  the  adjacent  electrode  made  cleaning  more  dif- 
ficult. The  acetate  cleaning  technique  cannot  be  extended  over 
the  central  area  (outside  of  the  re;gion  to  be  coated)  because  of 
the  vertical  extension  of  the  electrodes,  particularly  at  the 
attachment  points.  The  central  area  is  not  as  clean  as  the  sur- 
rounding area  and  contamination  may  result  if  particles  become 
detached.  Having  no  success  with  the  acetate  sheet  method  a 

1 1 

series  of  wash  and  soft-wipe  procedures  was  used. 

1 1 

Filtering  the  quantity  of  solution  required  for  coating  the 
large  disc  in  a short  time  (to  minimize  the  evaporation  of  the 
solvent  during  application)  created  problems.  The  high  pressure 
syringe  caused  some  solution  to  leak  from  the  Luer  lock  of  the  syringe 
filter  connection.  Other  lower  pressure  Teflon  filters  were  used 
but  these  had  produced  even  worse  leakage  problems  within  the 
filter  holder.  The  presence  of  a leak  did  not  necessarily  con- 
taminate the  coating  if  precautions  were  taken;  however,  control 
of  a leak  diverted  attention  from  the  coating  effort.  Recencly, 
new  nylon  filters  were  tried  on  smaller  coatings  and  problems 
with  leaks  and  high  back  pressure  were  corrected,  Figure  16. 

With  a smaller  coating,  the  intrusion  by  the  coating  appara- 
tus was  shorter  in  time,  and  less  extended  physically  requiring  that 
only  the  stir  rod  and  the  filter  apparatus  enter  into  the  still 
region  above  the  solution  layer.  With  the  larger  coating  the 
syringe,  the  hands,  and  arms  of  the  coater  were  directly  in  that 
region.  Care  was  taken  so  that  these  additions  do  not  introduce 
particle  matter.  Precautions  included  use  of  a filter  mask  over 
the  nose  and  mouth,  nylon  cap  and  a coat  with  disposable  gloves 
taped  to  the  sleeves  of  the  coat  , Figure  14..  The  fact  that  these 
parts  extended  over  the  coating  region  and  are  there  for  a 


FIGURE  16.  COATING  TOOLS  (The  precision  level,  upper  left,  a 
1 Op  nylon  filter,  lower  left.  Guth  wash  bottle,  upper  -ight 
used  for  cleaning  other  instruments.  Filter  holder  (with  attach 
syringe)  for  TP  solution  cleaning,  middle  right.  Stir  rod,  lower 
right. ) 
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longer  period  of  time  during  coating  does  increase  the  risk  of 
contamination . 

After  the  coating  procedure  is  completed  and  the  coating  has 
set,  it  is  placed  in  a vacuum  bell  jar  and  is  heated  to  about  50  C 
for  at  least  24  hours  to  drive  off  all  remaining  solvents. 

The  solvent  used  was  48  per  cent  toluene,  32  per  cent  methyl 
ethyl  ketone,  and  20  per  cent  ethyl  acetate.  The  thermoplastic 
was  then  dissolved  in  this  mix  in  an  amount  which  measures  six 
per  cent  by  volume.  Typically  6.6  cc's  of  solution  was  coated 

O 

onto  an  area  of  600  cm  to  produce  a solution  layer  110  microns  thick 
which  upon  evaporation  produces  a thermoplastic  layer  about  6%  mi- 
crons thick.  The  solution  was  filtered  with  a 1.0  micron  nylon 
filter,  type  "NR"  produced  by  Millipore  Corp. , Bedford,  Mass. 

An  example  of  coating  results  is  gi'  en  in  Figure  17  which  is 
a composite  (10  circular  areas)  of  Mach-Zehnder  interf erograms  of 
the  coated  15"  diameter  disc.  The  fringes  were  intentionally  set 
in  during  the  Mach  Zehnder  testing.  The  straightness  and  parallelism 
of  the  fringes  indicates  very  good  coating  and  substrate  quality. 
Thickness  uniformity  over  2"  regions  is  better  than  one-fourth 
wavelength. 

The  coating  method  developed  for  the  smaller  5"  diameter  sub- 
strate is  similar  to  that  discussed  above  and  has  been  described  in 
a previous  report.  Interf erograms  of  the  5"  substrate  are  included 
here  for  completeness  and  are  shown  in  Figure  18.  In  Figure  18a 
the  interf erogram  is  set  for  the  minimum  possible  number  of 
fringes  which  can  be  seen  to  be  zero.  To  allow  recognition  of 
fractional  wavelength  nonuniformities  not  observable  in  18a  fringes 
have  been  included  in  Figure  18b  which  is  a repeat  of  18a 
except  for  the  Mach-Zehnder  fringe  pattern  setting.  Here,  better 
than  sixth-wave  uniformity  has  been  achieved  over  a central  1" 
square  region. 


FIGURE  17 . COMPOSITE  INTERFEROGRAM  OF  A TP  COATED  15"  DISC 

(~one-fifth  size) 
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RECORDING  FORMAT 

The  mathematical  formulation  for  two  dimensional  Fourier  trans- 
formations of  phase  modulation  functions  (signals)  in  polar  and  in 
rectangular  coordinates  is  reviewed  herein.  The  case  of  sinusoidal 
modulation  which  is  synchronized  with  the  start  of  each  scan  is  treated 
in  Sections  6.1  and  6.2  for  rectangular  and  polar  formats.  Unsynchro- 
nized modulation  is  analyzed  in  Section  6.4.  Additional  analysis  of 
polar  format  recording  is  given  in  Sections  6.3  and  6 . 5 by  methods 
which  allow  a more  comprehensive  treatment  of  the  unusual  properties 
of  the  Fourier  transform  of  polar  data  not  fully  covered  in  Section  6.2 


r 


FIGURE  19.  COORDINATE  GEOMETRY 

The  coordinate  geometry  is  defined  in  Figure  19.  In  polar  coordinates 
the  input  signal  is  defined  in  the  r,  6 plane  and  is  denoted  as  g(r,6). 
The  Fourier  transform  of  g is  denoted  as  G(p , $)  and  is  defined  in  the 
p,  <j>  (frequency)  plane.  Similarly,  in  rectangular  coordinates,  the  in- 
put signal  to  be  transformed  is  g(x,y)  and  it  occurs  at  the  input  x,  y 
plane.  Its  Fourier  transform  G(wx,  ay)  is  in  the  frequency  plane  which 
has  the  space  variables  u,  v. 

The  input  data  g is  assumed  to  be  a light  field  which  is  realized 
by  coherent  optical  readout  of  recorded  data  which  has  been  generated 
as  a line  scan  recording.  In  the  polar  format  each  scan  line  path 
occurs  in  the  radial  direction  and  adjacent  scan  lines  are  spaced 
by  an  angle  6 with  each  line  having  a finite  width  0T.  The  re- 
corded signal  is  modulated  along  each  radial  scan  line.  A polar 
format  signal  with  the  modulation  cosu^r  along  each  scan  line  is 
sketched  in  Figure  20. 
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FIGURE  20.  INPUT  SIGNAL  FORMAT 
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Figure  20  depicts  a line  scan  recording  having  modulation 

cosoj^y  along  each  scan  line  in  the  rectangular  format.  In  this 

case  the  adjacent  scan  lines  are  spaced  a distance  and  each 

line  has  a width  x . 

x 

The  case  of  interest  here  is  that  of  a coherent  optical  Fourier 
transform  channel  where  the  input  data  exists  in  the  form  of  vari- 
ations in  the  complex  optical  transmission  t = TeJ  3 of  the  re- 
cording medium.  In  particular,  the  recorded  data  will  be  in  the 
form  of  a surface  deformation  modulation  denoted  by  deformation 
depth  d.  The  deformation  modulation  will  affect  the  phase  a of 
the  optical  transmission  function  t.  The  transmissivity  T is 
constant  (T  = Tq)  . Upon  readout  with  a collimated  coherent  light 
beam  which  illuminates  the  area  of  the  input  data  to  be  processed 
we  realize  for  our  input  signal  a light  beam  emergent  from  the 
recording  medium  which  is  equal  to  the  product  of  the  illumina- 
tion light  field  amplitude  E and  the  complex  optical  transmission 
t.  Thus,  the  input  g,  can  be  expressed  generally,  in  polar  form, 
as 


g(r,6) 


jr'j  gju[d(r,0)] 
o 
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and  in  rectangular  form  as 


g(x,y)  = ET0ej“td<x’3,>I 

We  will  be  taking  E as  constant  over  the  input  aperture. 

The  recorded  signal  deformations  d are,  of  course,  related  to 
the  internal  workings  of  the  recording  device.  Assuming  the  re- 
cording device  to  be  comprised  of  an  electron  beam  which  writes  on 
a deformable  thermoplastic  then  d would  be  proportional  to  the 
recorder  input  that  is  used  to  modulate  the  writing  electron  beam. 
The  performance  of  the  recording  device  will  not  be  considered 
explicitly  here. 

At  the  output,  or  frequency  plane,  we  obtain  the  optical 
Fourier  transform  G of  the  input  g,  i.e.  , G = F [g]  . The  transform 

in  polar  coordinates  for  polar  format  input  data  can  be  expressed 
as 


G(p,<p)  = A Idr  *r  Id0  *g(r , 0)  e"^7rrcos  ^ 0 “ ^ 


It  should  be  noted  that  this  expression  represents  a Fourier  trans- 
form in  the  p dimension  only  because  it  is  based  upon  optical  pro- 
cessing elements  having  circular  symmetry  such  that  no  optical 
power  occurs  in  the  circular  dimension. 

In  rectangular  coordinates  we  have  the  Fourier  transform  ex- 
pression 


dx  dy 
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This  rectangular  coordinate  expression  is  also  based  on  optical 
processing  elements  having  circular  symmetry  as  with  a spherical 
lens  for  example. 

The  term  A appearing  before  the  transform  integral  contains 
the  amplitude  constants  and  the  phase  terms  of  the  coherent  opti- 
cal Fourier  transformation  process  which  are  not  dependent  on  the 
variables  of  integration. 

6,1  RECTANGULAR  FORMAT  DATA 

The  input  g(x,y)  to  be  transformed  is  the  rectangular  format 
line  scan  phase  recording  with  sinusoidal  modulation  in  the  scan 
line  direction  y,  as  sketched  in  Figure  20.  The  scan  line  de- 
formation grooves  in  the  y direction  are  periodic  in  the  variable 
x with  period  x and  width  x„ . The  x direction  deformation 

p T 

variation  is  described  by  a function  h(x)  that  is  normalized  re- 
lative to  the  peak  deformation  Dx  and  can  be  expressed  generally 
as  a Fourier  series 


h (x)  = £ C.  cos  TT-  x 
k K x- 


2TTk 

t 

P 


where 


ck* 


-*J 


xt/2 


■x, 


h(x)  cos  I xjdx 

\ x 

,/2  ' P 


/2irk  \ 

Mc 


A sinusoidal  modulation  in  y along  each  scan  line  is  denoted  as 


f(y)  = D cos  w^y 


The  lWO  dimensional  deformation  expression  is  then 
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The 


d(x,y)  - (D  + D cos  u>-y)  £ n cos  kw  x 

1 k k 0 

where  is  the  scan  line  spacing  frequency,  i.e. , w = ^ 
input  signal  to  be  transformed  is  now  ° XP 


g(x,y) 


; 2 it  (n  - 1) 


ET  e 
o 


(D  + 
x 


D cosw 


cos  ku)  x 
o 


where 


a 


2ir(n  - 1) 

X 


D 


x 


b - jily  i)  D/2 

The  Fourier  transform  G(w  , w ) is  now  written  as 

x y 
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;i/z  r1'  E c 

G (ux>  wy)  - B f f e [ k 

-y^/ 2 -x^/2 


j fa  E cos  kcoQx  + b E c 


k cos  (ku)ox  + a)iy)J 


-j  (to  x + to  y) 

• e x y dx  dy 


y-,/2  x-,/2 

j*  ^ /*  -L 


/.!  /.I  tt  £+m  j £kto  x jm(kco  x±(o1y)n 

f I V R ('J)  J*(aVe  J.(bV*  ° 1 

-yi/2  -xx/2  R LA*m  J 


-j  (to  x + to  y) 

•e  y dx  dy 


If  we  take  the  scan  line  cross-section  h(x)  to  be  simply  a sinu- 
soid C Cos  toox  rather  than  a general  series  Eck  Cos  k to  x,  as 
is  very  nearly  the  case,  we  get 


_ X-.  sin  —ry  (to  + ( £+m)  to  ) 

G(tox,toy)  = B E J.9(aC)  Jm(bC)  4 2- 


2 x., 


~T  ^wx  + w0) 


yp 

y i sin  -j-  (w  ± mu,) 

• -ry  y 

^ yi 

~2~  (tOy  ± tom 


Here  we  have  used  a rectangular  input  data  aperture  over  -x^/2  < x 
x-^/2  and  -y^/2  £ y £ y 2 and  have  collected  non-integrand  terms 
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The  light  intensity  distribution  |G(u)x>  w ) | which  is  available 
for  observation  in  the  transform  plane  is  expressed  as 


9 9 ^1  ^1  ^ 

|G(o)x>  o)y)|2  = I B | 


^ J2(aC1)J2(bC1)Sinc2Ca)x+(^.4m)£a)o) 


£,m 


X1  2 i 

• j—  Sine  (to 
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for  cases  where  there  is  insignificant  overlap  of  orders.  Thus 
we  get  spots  of  light  in  the  frequency  plane  located  at  (tox>  w ) 
with  <ux  = (£+m)ku)o  and  oi  = mco^.  The  peak  intensity  is  propor- 
tional to  recorded  deformation  depth  through  the  arguments  a and 
b in  J [ ] and  Jm[  ].  The  spatial  distribution  of  the  light  spots 

is  given  by  the  Sine  function.  In  practical  use  for  spectrum 
analysis  of  the  recorded  data  we  choose  recorded  depths  to  retain 

mainly  two  pairs  of  spots  of  light  in  the  frequency  plane  (w  = 

+o)^,  wx  = wq)  for  any  one  recorded  input  frequency  co^. 

6.2  POLAR  FORMAT  DATA 

The  input  g(r,  9)  to  be  transformed  will  be  the  polar  format 
line  scan  recording  with  sinusoidal  modulation  in  the  line  direc- 
tion as  sketched  in  Figure  20.  The  scan  lines  are  straight  de- 
formation grooves  in  the  r direction  which  will  be  periodic  in 
the  polar  angle  variable  9 with  period  0^ , groove  width  9^  and 
mean  groove  depth  Dg.  The  angles  9,  0 , and  9^  are  in  radians 
and  Da  in  millimeters.  The  scan  line  is  of  constant  width  in 
linear  measure.  Its  angular  width  0^  will  vary  inversely  with 
r.  The  scan  line  description  in  the  0 direction  will  be  denoted 
by  a function  h which  is  normalized  relative  to  the  peak  depth 
Dq . Since  it  is  periodic  in  0 we  can  express  h with  the  following 
Fourier  series 


0 


h(r , 0)  = E C,  (r)  cos 


2nk 


where 


C,_  = 


k e 


7- 


h(6)  cos  0 d 0 


Tie 


Deformation  modulation  along  the  scan  line,  denoted  as  the  func- 
tion f will  be  sinusoidal  with  peak  deformation  D and  frequency 
, thus , 


f(r)  = D costo^r 


The  composite  recording  in  r and  0 is  written  as 


0 


d(r,0)  = (D0  + D cosw^r)  £ C^Cr)  C0S  ^""0 

k 


2tt 


where  tofl  is  the  scan  line  spacing  frequency,  i.e.,  = j—.  The 


signal  to  be  transformed  g(r,0)  is  expressed  as 


g(r , 0)  = ETQe  L 


2TT(n  - 1) 


A 


(D0  + D cosco  Xr)E  C^.(r)  cos  kco0( 

k 


= ET  e 

O L 


jja  ^ (^(r)  cos  kw00  + b ^ cos(ku)Q0+  oj^r) 
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D/2. 


where  a 


2 it  (n-1)  n 

X u0 


and  b 


2 tt  (n-1) 

X 


Using  the  Bessel  series  expansion  for  the  exponential  terms  in 
g(r,0),  its  transform  can  be  expressed  as 


G(p,4>) 


- * 1+m+n 

“ Bjdr  ryd0-  J t ^ VaCk)Jm(bCk)Jn(cr)e 


jx-  .R  0)^0 


k X,m,n 


jm(kaj00  ±3r)  -jn(0-<J>) 


where  B = AET  , c = 2 tt p = w.  and  the  limits  of  integration  for 
0 0 0i  0 ri 

0 are  0 - ^ to  0 + -^  and  for  r are  r - to  r + . 

o ~2  2 o 2 o TT 

We  will  continue  with  the  assumption  that  the  scan  line 

cross  section  is  sinusoidal  in  0 giving  C Cos  ui00  rather  than 

C^  Cos  kwQ0.  The  integration  on  0 is  then  as  follows. 


j [(X  + m) 
e 


u -n]  0 jn<|> 

0 dee 


-j  (X  + m) 
e 


we 


1 

T 


. 


where  q = ( £+m)  w 


0 


n 


1 2tt 

We  must  note  here  that  to0  = 2tt  — where  0 = — - and  N is  very 

4 P ^ 

large  ( < 10  ) . The  above  result  contains  a Fourier  series  expan- 
sion over  q which  defines  a pulse  of  width  0^  in  the  <j>  dimension  of 
the  frequency  plane.  The  pulse  has  period  2 tt  and  can  have  fine 
structure  within  the  pulse  depending  upon  the  Sine  function  be- 
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havior  over  the  index  variable  q.  The  pulse  in  <j>  has  width  3 
which  is  identical  to  the  input  plane  pulse  width  6 . The  polar  for 
mat  Fourier  transformation  now  has  the  form 


Be 


G (p  , 40  _ 2 


x 

i fr  +' 


dr  *r* 


E 

i ,m,n 


, . N£.-fm+nT 
(-j)  J 


£^aC  9)  Jm('bC  y)  jn('cr^ 


jmw-^r  — jncj> 
• e e 


sir.  q 
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We  will  continue  with  the  integration  on  r with  the  approximations 
that  Sine  (q)  is  separable  from  the  series  on  and  that  the 
integrand  term  r has  a minor  effect  and  is  taken  constant  over 
the  range  of  the  integration.  With  these  assumptions,  we  can 
write  G as 

D <-j>i*v*c>  Vbc>  s-^r 

£,m,q  q-  2 


+jmw,  r 

dr  El  (-j)nJn(cr)e 
n 


The  assumption  that  Sine  (q)  is  separable  from  the  Jn  series  appears 
reasonable  because  the  value  of  u)Q  is  very  large.  With  the  separability 
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assumption  the  Jn  series  exists  alone  under  the  integral  operator 
and  can  be  changed  as  follows: 


!(-j)nJn(cr)  = e^cr,  where  c = 2np  = to 


The  integral  on  r has  the  form 


/ 


ro  + 2 ej(c±mus)rdr 


sin  (a L±ma),  )-T 

ej(a3p±muil)ro 


ro  ' ~T 


(0)  im^)^ 


This  approximation  for  the  Fourier  transform  G(p,<J>)  is  thus 

6, 


Beiri 


G (p  , 4»)  - — ^ 


X)  (-j)£+mJ£(aC  )Jm(bC  )S1!qT"  e"jn<J> 
£,m,n 


01 

qT* 


rx 

sin  (to  ±mto,  )-*-  s 

. _e ej(uip±inal)ro 

rl 

(to  imo^)-^ 


B0lrl 
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ri 

_ ~ sin(  to  +mto,  )-rr  , v 

S g (4>)  Rect  ($/ 0^)  £ — eJ(wp-mwi)  o 

m ( to  ±mto^)-7jr 


We  note  that  the  combined  terms  J^J^e-^  a)0~q^sin  q-^  v -j1  q 

are  a Fourier  series  for  the  original  g(0)  • (with  costo-^r  = 1) 

o 

but  now  in  the  variable  tj>,  i.e.,  g(<f>)  taken  over  the  interval 
0 1 0 1 

- y < <f>  < ~y  as  indicated  by  the  rect  (<J>/0-^)  function. 


The  light  intensity  distribution  |G(p,cJ>)|  is  available  for 
observation  in  the  transform  plane  and  it  is  expressed  as 

(0  r \ 2 

Rect(<f)/ep  Y sine2 

for  cases  where  there  is  insignificant,  overlap  of  orders.  We  get 
arcs  of  light  with  radial  positions  u.  = ± mu and  the  arcs  have 
angular  width  0-j^  given  by  the  rect  function.  This  approximate 
solution  is  an  incomplete  representation  of  what  can  be  demon- 
strated experimentally.  The  main  difference  is  the  absence  of 
terms  representing  higher  order  diffraction  in  the  4>  direction. 

An  alternative  approach  to  the  analysis  of  polar  format  data  is 
given  in  Sections  6.3  and  6.5. 

6.3  EXPERIMENTAL  ANALYSIS  OF  POLAR  DATA 

It  is  important  to  review  experimental  results  for  the  optical 
Fourier  transform  of  polar  formatted  data  as  a guide  when  deriving 
a mathematical  representation  of  frequency  plane  distributions. 
Experimental  analysis  of  data  recorded  in  a polar  format  was 
accomplished  with  inputs  of  the  form  shown  in  Figure  2l  which 
consists  of  two  adjacent  30°  sectors  each  with  synchronized  pulse 
modulation.  The  ratio  of  the  pulse  repetition  rate  of  sections 
A to  B is  5 to  3 . Other  specifics  are  given  in  the.  figure. 

The  first  few  transforms  to  be  considered  are  apertured 
sections  of  the  transparency.  The  aperture  was  a 7.5°  angular 
sector  as  illustrated  in  Figure  21b,  When  placed  over  section 
A,  the  transform  shown  in  Figure  22  was  obtained,  and  when  over 
i B,  the  transform  shown  in  Figure  23  was  obtained.  The  vertical 

direction  f^  corresponds  to  the  centerline  direction  y of  the  input 
aperture  (see  Figure  21b).  A study  of  these  transforms 
I 

I 

1 

1 
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FIGURE  22 . FREQUENCY  PLANE  FOR  INPUT  APERTURE  A 


indicates  that  pulse  modulation  in  the  input  plane  produces 
characteristics  in  the  transform  plane  that  are  discussed  below 
and  shown  graphically  in  Figure  24 

1.  A multiplicity  of  orders  for  which  each  "center  of  mass" 
is  on  a rectangular  grid  with  coordinates 


( 


M 

6 R 

p mean 


where 

M & N are  the  order  number  as  indicated  in  Figure  24 
Rb  + R 

*\nean  ~ 2~  ■'■s  t^e  mean  position  of  the  aperture, 

Rp  is  the  nulse  modulation  period, 

0p  is  the  azimuthal  period  of  the  scan  lines 

Each  order  is  confined  to  an  area  about  its  center  of 
mass,  where  the  area  is  distorted  with  respect  to  the 
input  aperture  and  where  the  distortion  depends  on  the 
order  number. 


2. 


3. 


The  orders  along  the  f direction  are  contained  in  angle 
eA  about  fy»  where  0^  is  the  input  aperture  angle 
These  orders  are  the  FT  of  the  pulse  modulation,  e.g. 

the  orders  occur  for  f = — ^'■e&er  f and  their  width  along 

y P 

f is  inversely  proportional  to  the  scan  line  length  R,  -R 
y b a ' 

The  orthogonal  direction  f^  represents  the  frequency  con- 
tent or  shape  of  the  scan  lines.  If  we  define  a tangen- 
tial frequency  w = in  the  input  plane  (i.e.  the  spatial 
frequency  of  the  scaR  lines  along  a circle  of  radius  R) 


then  the  Mth  order  (along  fx  where  N = 0)  is  confined  to 
the  interval  - to  — . That  is,  it  is  confined  to 


0 R ‘ 
p a 
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the  extremum  tangential  frequencies  of  the  input  and  there 
is  an  inverse  '•elationship  of  the  nature  that  if  the  aper- 
ture Is  closed  down  from  the  outside  (i.e.  reduced) 
then  the  orders  along  f reduce  in  size  from  the  inside. 

This  means  that  the  orders  along  f will  eventually  over- 

X 

Ilap  since  their  size  increases  with  the  order. 

We  also  note  that  these  orthogonal  orders  are  con- 
tained within  the  aperture  angle  0^  centered  along  fx  axis 

4.  There  does  not  seem  to  be  any  well  defined  fringe  structure 
within  any  order  except  when  orders  overlap.  This  latter 
result  is  seen  in  Figure  25  for  which  the  aperture  of  Figure 
21b  was  centered  on  the  boundary  between  both  sections,  i.e., 
3.25°  of  section  A and  3.25°  of  section  B.  Regions  of 
overlap  between  the  two  sets  of  orders  for  the  two  differ- 
ent sections  produce  fringes  indicative  of  the  pulse  modu- 
lation ratios  of  the  two  sections.  Figure  26  is  for  the 
same  kind  of  setup  as  for  Figure  25  but  where  the  radius 
R of  the  aperture  was  reduced  to  2 mm.  This  figure  shows 

Si 

both  the  interference  effects  between  the  different  orders 
and  the  inverse  characteristic  of  the  orthogonal  orders 
described  in  (3)  above.  The  flaring  seen  at  the  end  of 
the  orders  is  due  to  the  smearing  together  of  the  scan 
line  and  pulse  modulation  at  the  smaller  radii  on  the  trans- 
parency. 

Finally  the  transform  for  the  entire  transparency  (60°  sector) 
is  shown  in  Figure  27  . There  is  significant  overlap  of  many 
orders.  If  the  FT  orders  along  the  f^-direction , within  the  0-j. 
angle,  are  to  be  clear  of  any  other  orders  then  from  geometrical 
consideration  in  Figure  24  we  need  to  specify  that 


or 


arctan(0pVmax) 


-1 


where  r is  the  maximum  spatial  frequency  of  interest, 
max 

An  alternative  to  the  above  criteria  may  be  to  eliminate  the  higher 
order  M,  N > 1 by  using  a scan  line  with  a sinusoidal  cross  sec- 
tion that  linearly  expands  with  R such  that  there  is  only  a single 
azimuthal  frequency.  It  has  not  been  shown  experimentally  that  such 
a condition  does  indeed  eliminate  the  higher  orders. 

It  should  be  pointed  out  that  many  of  the  observa- 
tions made  are  dependent  on  the  input  aperture  being  not  too  small 
(in  order  to  avoid  significant  aperture  diffraction  effects)  and 
not  too  large  (in  order  to  avoid  overlapping  orders). 

6.4  UNSYNCHRONIZED  MODULATION  COUPLING 

The  example  of  sinusoidal  phase  modulation  along  scan  lines, 
in  either  polar  or  rectangular  format,  was  described  in  the  pre- 
vious sections  in  a way  which  assumed  that  the  phase  of  the  sinu- 
soid was  the  same  at  the  start  of  each  scan,  i.e.,  synchronized. 

When  this  condition  is  not  satisfied  then  frequencies  can  be 
generated  in  the  direction  orthogonal  to  the  scan  lines. 

In  the  rectangular  format,  if  we  consider  some  arbitrary 

modulation  frequency  f^  in  the  scan  line  direction  and  a scan 

line  repetition  frequency  f , then  it  can  be  shown  that  we  can 

s 

realize  synchronized  recording  if  f^  is  an  integer  multiple  of 
f . That  is,  if  q is  an  integer,  we  require  for  synchronization 


when  the  above  relationship  does  not  hold  true  and  is  some 


value  between  (q  + l)f  and  qf  we  can  define  the  value  of  fre- 

s s 

quency  induced  in  the  direction  orthogonal  to  the  scan  line  as 
follows.  First  we  examine  f^  relative  to  fg  and  determine  num- 


bers n and  a such  that 


q - <n  + a>  fs  * nfs  + afs 


with  the  constraints  that  n = positive  integer,  1 < a ± 0. 

Then  the  resultant  spatial  frequency  component  f in  the  di- 
rection orthogonal  to  the  scan  lines  (x  direction  for  the  rectangu- 
lar format)  is 


af  y 
f = -511 
a x_ 


Here  f, , f and  f are  spatial  frequencies,  y is  the  spatial 
las  P 

equivalent  of  the  temporal  scan  period  and  x^  is  the  spatial 
spacing  between  successive  scan  lines. 


A similar  result  exists  for  the  polar  format  except  that  we 
work  in  the  variables  r and  9.  Thus,  in  the  polar  format  we  have 


af  R 


fa  * 


We  must  introduce  this  added  modulation  to  the  input  signal 
functions  for  the  respective  cases  of  g(x,y)  and  g(r,0).  It 
must  be  included  as  a phase  function  which  is  added  to  the  sinu- 
soid of  frequency  f^  which  was  considered  in  the  above  discussion. 


, ft  . ■ 


V ,:y.M L.  ; . . 


to— -a 


Thus  in  rectangular  coordinates,  the  scan  line  modulation 
D cos goes  to  D cos  (to^y  + ^x)  . In  polar  coordinates  D cos 
(J^r  goes  to  D cos  (oo^r  4-  to  0)  . 

Regarding  resolution  in  the  frequency  plane  we  note  (in  the 
rectangular  format  variables)  that  a spot  of  light  in  the  fre- 
quency plane  will  be  located  along  loci  which  are  a set  of  parallel 
sloped  lines  in  the  to  , cox  plane.  This  can  be  seen  by  examining 

the  Fourier  transform  expression  G(oo  , to  ) when  the  to  term  is 

x y s. 

included.  The  slope  of  the  loci  in  the  <o  , oox  plane  is 


x 

-2-  = slope  of  loci 


The  set  of  loci  cross  the  to^  axis  at  the  frequency 


CO 

y 


q(jjg,  loci  crossing  on  to  axis 


The  line  spacing  of  the  loci  in  the  co^  and  tox  directions  are 


co  = loci  spacing  in  to  direction 

s y 


ws  SE*  = loci 
P 


spacing  in  to  direction 

X 


The  width  of  a single  frequency  spot  in  to  and  to  is  defined  by 

y x 

the  sine  function  expressions  of  the  Fourier  transform  G(to  , to  ) 

x*  y 


55 


k 


where  we  find  the  spot  half  width  from  its  center  to  its  first 
zero  is 


Aw  = — - . for  w direction 

y Y]_  y 


Aw  = — , for  w direction 
X X X 


And,  for  comparison  the  ratio  of  spot  width  to  loci  separation 
in  the  two  dimensions  is  seen  to  be 


w y1 
sJ  1 


= 1,  for  Wy  direction 


w y x, 
s^p  1 


= — , for  w direction 


Improved  resolutions  in  the  w^  direction  is  desirable  and  may 
be  available  by  use  of  double  raster  recordings  . 


6.5  FURTHER  DERIVATION  OF  POLAR  DATA  TRANSFORM 

In  this  section  an  expression  for  the  Fourier  Transform  (FT) 
of  polar  formatted  data  will  be  discussed  which  describes  higher 
order  diffraction  in  both  the  dimensions  of  the  frequency  plane 
p,  <f)  rather  than  only  one  dimension  as  described  in  Section  6.2. 
This  analysis  assumes  recorded  data  which  can  be  represented  as 
optical  transmission  variations.  This  will  be  true  for  phase 
modulation  where  recorded  deformation  depths  are  small. 


. 
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First  we  will  consider  a series  of  polar  scan  lines  with 
no  modulation  along  the  scan  line  in  order  to  determine  the  FT 
of  the  basic  system.  A single  scan  line  of  length  = AR 

and  constant  width  W is  described  by  (see  Fig  28^ 


recti— Jrect 

W 


(ir) 


*<$(y  -ARC) 


R + R. 
where  R = — -a 


Then  its  FT  is  WAR  sir.c(Wfx) sine (ARfy)e~^Tr'*"^y^'c 


The  transform  in  polar  coordinates  (p,(j>)  is 


WAR  sinc(Wpcosc}>)sinc(ARpsincf>)e 


-2TriR  psintj) 


If  there  are  N scan  lines,  separated  by  the  azimuth  period  0 , 
then  the  total  signal  is  (see  Fig  28b) 


G(p,cj>)  = WAR  XI  sine  [wpcos  (c|>+n6  )] 
N-l  p 


sine  [ ARpsin  (4>+n0p)  ] 


-2iriRcpsin  (cj)+n0p) 


where  for  convenience  we  let  N be  an  odd  integer.  This  equation 
(1)  is  our  basic  relationship.  Let  us  now  consider  the  condi- 
tions under  which  this  expression  is  to  apply.  First  we  have  that 
(see  Fig.  28b) 


e = Ne  < f 

A p 4 


mmm 


(2a) 


e <<  eA 

P A 


W <<  AR 


AR  ~ R 


However,  for  the  moment,  we  will  assume  that  both  W and  AR 
are  very  small  so  that  both  sine  functions  in  ( 1 ) are  nec\rly 
unity  and  can  be  removed  from  the  summation.  The  expression  we 
are  initially  interested  in  is  (see  Fig  29  ) 


-iMsin(^+nG  ) 
I<P,*)«Le  P 

n 


M = 2ttR  p 


That  is,  we  are  considering  the  summations  of  the  phase  terms  from 
a series  of  point  sources  at  a radius  Rc  separated  by  the  angle 


■i ng^ 


sinc^y^j  = 0 for  n f integer  times  N 


(5e) 


N-l 

~T~ 

sinc^(t)  ^ ' 6(f-n)  (N  odd)  (5f) 


The  full  width  of  sinc^t  about 
(to  the  first  zero) 


is  main  maxima 


18  = N 


(5g) 


All  secondary  sidelobes  are  less  than  ^ 


(5h) 


and  where 


00 

nr(t)  = £ 6(t~n) 
n=-°° 


( 5 i ) 


2ir 


We  now  let  0 = _ 

P p 

teger.  Then  (4)  becomes 


where  P is  a 


(large)  positive  and  odd  in- 


(6) 


k - mP  m an  integer 


which  defines  the  mth  diffracted  order, 
around 


That  is,  only  a few  terms 


J p(M)e-imP 


are  not  significantly  suppressed.  The  peak  value  for  this  term 
lLgl  fr01"  C°nSlderlng  the  characteristics  for  k 


maxJk(M)  7c  . 7k“1/3 


when  M as  k k >>  1 


half  width  of  this  main  peak  « k3^3 
area  under  main  peak  « 1.4 


(9c) 


Thus  taking  mP  . M to  satisfy  condition  (9a)  we  have  from  eqs 
( ) and  ( b ) that  the  Bessel  function  of  order  mP  peaks  at 

pn  and  has  a width  Ap  where 


p = _SL_ 

m 0 R 


' 


» V - Mj  • •.  . i : 


I 


I 


l 


Ap  = 


Pm 


m 


(mP) 


YTS 


(ID 


m 


is  m 


times  the  tangential  frequency  of  the  point  sources  and 


Ap  is  a small  portion  of  p since  P is  a large  number, 
m r m 


th 


Let  us  now  find  the  extent  of  the  mLa  order.  Recall  from 
(2a  ) that  N < P.  We  let  P/N  = Q where  we  assume  that  0 is  an  odd 
integer.  Note  that  0.  = 2tt/Q.  Hence,  as  indicated  in  Figure  30, 
there  are  Q terms  contained  within  the  half-width  of  the  main 


maxima  of  the  sinc^  function. 


t"H 

Thus  for  the  m order  the  signal 


is 


mP+ 


¥ 


E 


.7(mP)"1/3elk<)) 


(12) 


k = mP- 


0-1 


= .7(mP)“1/30  elmPc|)sinc 


A. 

0 2 TT 


Since  Q = then  the  sinc„  function  limits  the  mLa  order  energy 


th 


to  within  t^e  angles  <j)  = 4-  Note  that  due  to  the  periodicity  of 

2ir  for  the  sincQ  function  that  all  its  primary  maxima  fall  on 


top  of  each  other.  From  general  symmetric  considerations  for  FT's 
we  know  that  a symmetrically  placed  order  to  m must  exist.  It 
comes  from  placing  the  negative  value  of  m in  eq  (7).  From  eq  (12) 
we  also  see  that  the  orders  slowly  die  in  intensity  by  m ^3. 


The  range  of  p can  be  found  from  the  ranges  of  indices  allowed 
by  sincN  since  the  Bessel  functions  main  peak  width  greatly  ex 


ceeds  that  of  sincN  (see  eq.  9 ).  The  index  range  is  from  (mP 


to  (mP  + ^y^) , or 


- fei) 
2 ; 


i -J 


j 


1 

i 


i 


i 

i 


a 


64 


■ 


l 


m 1 ' n * m 4-  1 

^ - p - VP  ^Fc 


(13) 


h 

Thus  p is  limited  to  the  small  range  of  about  the  order  center 
which  implies  that  substituting  the  peak  value  for  the  Bessel  func- 
tion, as  we  did  in  eq  ( 12  ) , gives  quite  accurate  results  (see 
Figure  30)  . 

Now  at  this  time  we  can  go  back  and  consider  a set  of  N narrow 
scan  lines.  We  note  that  the  scan  lines  can  be  formed  from  the 


expression  ( 3)  by  integrating  over  R , that  is 


/ 


^b  -2friR'  psin(cj)+n0  ) 


P dR’ 


R 


(14) 


= A Re 


-2'rriRcpsin((j)+neT^) 


P sine  ^ARp sin  ((j>+n0p)J 


However,  rather  than  integrate  immediately  the  expression  ( 14)  we 
will  delay  the  integration  until  after  the  summation.  That  is,  we 
take  expression  (12)  to  be  approximated  by 


m 


c(4>)  5 (p - q ri1 


P c 


1 1 


l 


i ] 


1 


1 

! 1 

I 

I 


1 


I 


^ ZZ 


f 


If  AR  <<  Rc  then  pffi  can  be  substituted  for  Pm> 

Thus,  the  effect  of  the  sine  term  in  (14)  is  to  disperse 
the  spectrum  orders  proportional  to  the  order  number.  Hence 
the  orders  will  tend  to  overlap  at  high  order  numbers.  We  would 
expect  the  first  sine  term  in  (1)  to  also  disperse  the  signal, 
but  since  W is  small  its  effect  is  not  expected  to  be  significant, 
except  to  reduce  the  intensity  of  the  higher  orders. 

We  would  now  like  to  go  back  to  eq.  ( 3c  ) and  allow  a sinu- 
soidal modulation  of  the  point  sources  such  that 


N -iMsin(<j)+n0  ) 

I (p  , <t>)  = £ cos(en)e  p 

n 


= Jk<M)[sincN(l  + tf)  + sincIj  - 


Thus,  using  the  same  techniques  as  before  we  find  that  the  m 
order  is  now  split  and  given  by 


th 


"m+  9 R (m  ± 2tt) 

me'  / 


We  note  that  the  order  will  move  by  one  whole  order  position  for 
e = +2tt  . 


Modulated  Scan  Lines 

Let  us  now  consider  a sinusoidal  amplitude  modulation  of  the 
scan  lines  of  the  form  1 + Acos(2tTFr  + ne)  where  F is  the  spatial 


frequency  and  e is  the  synchronization  factor.  It  is  straightfor- 


identical  except  that  the  second  sine  term  is  replaced  by 


sine  |ARp sin (<j>+n0 


+ j | sine  |ARpsin(<j>+ne  ) - aRfJ 


ine+27riR  F 


(16) 


+ sine  |ARpsin(<f>+n0  )j  + aRF  e 


-ine-2TriR  FI 


This  form  of  the  equation  is  convenient  for  computer  analysis  but 
not  for  analytical  understanding.  Rather,  we  start  with  eq.  (4  ) 
but  where  the  ^(M)  term  is  replaced  by 

cos  ( 2-rrFr+ne  ) Jk  ( 2tt p r ) dr 
a 


i r 

= zJ  1 


r - R 


rect 


AR 


£ Jk(2wpr)e2’lirFdr 


elne+ 


c.  c. 


(17) 


1 f -2iriFR  I 

= j |ARsinc(ARF)e  CJ  * Ak(F,p)eine+  c.c. 

where  ijj,  represents  convolution  with  respect  to  F and  A,  (F,p)  is 
the  FT  of  the  Bessel  function,  i.e. 


foo  # ^ 

J,  (2irpr)e  ^irF^ 


Here  T^(*)  is  the  k order  Chebyshev  polynomial  which  has  the 
following  properties  of  interest 


T^(cosa)  = cos  (ka) 


Tk(-x)  = (-:..)kTk(x) 


(19) 


T (x)  = 1 
o 


For  the  moment,  we  will  assume  that  AR  approaches  infinity. 
That  is,  we  substitute 


6(F)  = lim  ARsinc(ARF) 

AR^°° 


(20) 


for  the  sine  function  in  (17).  Notice  that  this  implies  that  the 
scan  lines  pass  through  the  origin  and  extend  to  +°°  • We  can 
then  write  for  (17) 


7T 

= Aq(F, p)cos (ka)cos (ne+jk) 


where  we  have  let 


(21) 


cosa  = — 

P 


(22) 


which  can  be  done  since  in  eq . (18)  the  rect  function  implies 
I P I (assumed  F is  positive).  We  note  from  (18)  that 


A (F,p)  = (F/2p)  J.(2rpr) 

o i~2 j o 

ttVp  - F 


Eq.  ( 6 ) now  becomes 


co  N 

I'  (P,  4>)  = NA  (F,p)  E E 

0 k n 


L 

M • 

f 


- ik  4>+n  0 ) 
• cos  (ka)cos  (ne+^k)  e P 


where  it  is  understood  that  all  summations  are  symmetric,  i.e 


CO  CO 

£ = £ 

’ n n=-°° 


We  now  use  the  identity  (see  Figure  31) 


GO  IV 

^ £ e-2TTik(A+nB) 


k n 


= £e  2TtikAsincM(kB) 


- £ £ 6 (A+nB-k) 
k n 


-7  -~-~T-.v 


to  write  (24)  as 


I'(p.<f0 


A 

o 


* 

<t> 


6(<j>+a+|-)  4-  6((j)+a- 


j)  + 6(<t>-a~7jr) 


(27) 


In  order  to  consider  the  physical  meaning  of  this  expression  we 
first  note  that  the  delta  functions  represent  a locus  given  by 


<t>  + a + j ~ 0 


(28) 


which  gi-^es 

j 

sintj)  = +cosa  (29) 

- 


or 


fy  = IF  (30) 


since  psin<j>  f and  pcosa  - F.  That  is,  these  are  straight  lines 
with  y-axis  intercept  +F.  When  convolved  with  the  rect  function 


V 


(with  respect  to  the  angle  $)  we  see  that  it  restricts  4>  to 
angle  0^  about  the  directions  (+a  + , which  of  course  varies 

with  p.  For  p = F,  that  is,  on  the  circle  of  radius  F,  we  have 
that  a = 0 and  hence  $ is  confined  to  0A  about  + \ . The  ttt 
function  represents  a high  frequency  fringe  system  whose  loci 
are  given  by 


<}>  + a,  + 2 _ 


= 0 , +1 , 


F = cos (m0  )f  + sin(m0  )f 
p x p y 


This  represents  straight  lines  tangent  to  a circle  of  radius  F 

and  with  slope  m9  . Similarly  the  loci  of  constant  phase  as 

represented  by  the  exponent  term  are  straight  lines  tangent  to  the 

same  circle  but  with  slope  m8p/e.  The  term  Aq(F,p)  limits  the 

energy  to  outside  the  circle  p = F and  gives  a signal  die  off  pro- 
2 2-1/2 

portional  to  (p  - F ) 

The  above  results  have  a s trai[ ntf orward  physical  interpreta- 
tion. A single  sinusoidal  modulated  scan  line  of  infinite  length 
has  a FT  that  is  composed  of  3 parallel  lines,  one  through  the 
origin  and  the  other  two  displaced  by  +F  from  the  origin  such  that 
they  are  tangent  to  the  circle  p = F.  Figures  32a  to  32c  show  various 
relationships  between  the  parameters  of  the  transform  lines  (T-lines) 
We  note  that  for  any  circle  with  radius  p >_  F the  T-lines  cut  the 
circumference  over  a range  of  6^  into  equal  increments  of  arc  length 
p9  . This,  of  course,  also  holds  true  for  F = 0,  which  represents 

no  modulation  on  the  scan  lines . If  we  now  consider  the  implica- 

til 

tions  of  eq.  (10)  for  this  latter  case  we  see  that  the  m order 


FIGURE  32 


occurs  where  2ito  R p = 2irm.  That  is.  the  m order  occurs 
P ^ 

along  a circle  of  radius  p where  the  arc  increments  between 
° m 

the  T-lines  (which  in  this  case  are  lines  of  constant  phase) 
are  2it  phase  differences  of  integer  order  m.  Similarly,  for 
F f 0,  we  will  see  that  the  higher  orders  occur  along  those 
circles  for  which  the  separation  criteria  holds,  i.e.,  along 

V/  2 9 

F + p^.  We  further  note  that  since  the 
phase  shift  is  not  really  linear  with  respect  to  n but  goes  as 
sin  (4)  + n8p)  , and  as  indicated  by  eqs.  (6)  and  (12),  only  Q 
terms  of  the  N terms  contribute  significantly  to  building  up 
the  energy  of  an  order,  while  the  other  terms  are  so  non-cynchro- 
nized  that  they  contribute  little  to  the  formation  of  an  order. 

Let  us  now  consider  this  in  some  detail.  Carrying  out  the 
convolution  with  respect  to  F for  the  general  form  of  (27)  gives 


r tt  “1  -2iriRcg((J,+i) 

I = sine  I ARg(<j)+  e * B(cj),e) 


T tt  1 -2iriR  Cg(c)5-J) 

sinc[ARg(ct,-J)J  e * B(cJ)7e) 


where 


g (4) ) = -27rsin^~ 


N . 

B4  , e)  = £ eine6  (4>-n9  ) 

n 


n / \ 

"(f)rect(4) e 


Obviously  the  maxima  for  the  sine  functions  occur  when  condition 
(28)  is  satisfied,  that  is,  these  maxima  occur  along  the  lines 

fy  ~ -F"  Rcp  > ARp  >>  1 and  a » 0 the  four  maxima  are  physi- 

cally separated  (this  does  not  occur  in  the  zero  order,  which  we 
will  later  consider  in  detail).  Considering  only  the  maximum  for 
which  <p  w ^2  - a we  can  make  the  approximation 


gC^-Tj-)  « -2psinu  sin (<j> ' / 2) 


where 


' = <p+a-ir/2 


The  term  of  interest  then  becomes 

N /4> ' -n6  \ 

22  Sine  2pARsina  sin^ 2.)  exp  t 


4iTiRcpsina  sin 


'<p ' -ne 


2 ^ )+  ne  (38) 


Comparison  of  ( 38)  for  e = 0 with  the  original  expression  (1) 
shows  that  it  has  exactly  the  same  form.  Thus  the  same  type  of 
analysis  as  for  the  original  section  will  hold  if  we  substitute 


P + P sina  = tip  - f 


<j>  + 4>'  = <j>  - a - tt / 2 


th 

Then  the  m order  radius  (for  AR  + Q)  is  given  by 


77 


f>n,-VF  - 


i - e / ■:.  tt\  2 

6 R ) 

pc  / 


and  for  AP.  <<  R the  order  has  a width  of  ARpm/Rc  along  the 
variable  psina.C  The  angle  4, ’ is  confined  to  ±0A/2  These  effects 
are  illustrated  in  Figure  24.. 

The  zero-order  has  to  be  considered  in  detail  by  itself  since 
it  need  not  satisfy  some  of  the  assumed  approximations.  In  fact, 
its  character  is  quite  different  than  that  of  the  higher  orders. 

The  energy  to  form  the  higher  orders  comes  from  the  diffraction  of 
the  grating  formed  by  the  scan  lines,  whereas,  the  zero-order  does 
not  use  this  diffraction  effect  but  rather  the  diffraction  from 
the  radial  modulation.  Consider  that  the  tangential  direction 
along  the  circle  p = F (i.e.,  the  first  order  position)  is  a direc- 
tion orthogonal  to  the  phase  variation  direction  (see  Figure  32), 
and  hence  the  phases  can  not  add  up  as  for  the  higher  orders . This 
can  be  seen  mathematically  in  that  the  higher  orders  are  due  to 
the  primary  maxima  of  the  Bessel  function,  whereas,  the  zero-order 
is  a resonance  condition  with  the  asymptotic  oscillation  of  the 
Bessel  function.  That  is,  we  substitute  the  asymptotic  approximation 


J (2-rTpr)  = —~zz  cos(2r,pr  - ^ - jj) 
K TT/pr 


into  this  general  form  of  eq.  (17)  to  get 


N 00  1.1c  ( 6 “Hn  0 ) v / r\  'n , 77  x 

I V V e P / k dr  cos  (2TrFr+nt.)  cos  (2irpr  -^c-^) 


nVpRc  n k 


where  we  have  let  l//r  = 1//RC  since  l//r  does  not  change  signi- 
ficantly over  the  range  R&  < r < R^.  Then  we  find  that 


.. 


T _ nar  . ATW„ 

I = smc  AR(F  - p) 

2tt/^R"  r 

c 


rect 


ft) 


\ 19  J 27ri(F-p)R 


+ c . c . * 6(4> 


Thus  the  zero  order  is  confined  to  ±-£  about  <j>  = + l and  has  a 
radial  width  of  about  i about  p = F.  We  note  that  the  peak 
intensity  at  p = F decreases  as  1/F  (while  the  arc  length  of  the 
order  increases  proportional  to  F keeping  the  total  energy  in  the 
order  constant) . 

Thus,  this  completes  the  construction  of  Figure  24,  which 
explains  the  experimental  results  given  in  Section  6.3. 
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6.6  OPTICAL  RECTIFICATION 

Analysis  of  the  spatial  frequency  spectrum  which  results  when 
data  is  recorded  in  a polar  rather  than  rectangular  format  indicates 
that  the  frequency  plane  distribution  may  not  be  suitable  for  a spec- 
trum analyzer  application  without  some  modification  in  the  processing 
optics.  Recall  from  Sections  6.3  and  6.4  that  a single  frequency  re- 
corded in  the  polar  format  results  in  a Fourier  frequency  plane  distri- 
bution that  is  spread  over  an  arc  rather  than  the  point-like  spot  ob- 
tained with  the  rectangular  input  data  format.  It  appears  that  fre- 
quency resolution  for  un-synchronized  sinusoidal  data  (Section  6.3) 
may  be  compromised  because  of  the  arc  distribution. 

Several  approaches  to  a correction  for  the  polar  format  may  be 
possible.  One  approach  considered  during  this  program  was  to  pro- 
vide, as  part  of  the  data  processing  optics,  optical  elements  which 
convert  the  polar  data  into  rectangular  data  (optical  rectification) . 

It  must  be  noted,  however,  that  since  coherent  optical  processing  of 
phase  data  is  being  considered,  the  light  field  amplitude  and  phase 
rather  than  its  intensity  must  be  properly  preserved.  This  is  in  con- 
trast to  some  extent  to  the  design  of  an  optical  rectification  scheme 
for  image  projection  devices. 

The  use  of  a spherical  lens  tilted  relative  to  the  recorded  data 
plane  allows  the  convergent  scan  lines  of  a polar  recording  to  become 
parallel.  However,  this  causes  the  ends  of  the  scan  lines  to  line-up 
on  a curve  rather  than  a straight  line  as  with  a rectangular  raster 
and  some  loss  in  resolution  also  occurs.  Added  provisions  for 
straightening  the  raster  edge  and  correcting  for  resolution  loss  were 
considered  but  the  methods  pursued  during  the  limited  effort  expended 
on  this  task  did  not  result  in  a satisfactory  solution. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Evaluation  of  the  performance  properties  of  the  experiment- 
al model  real-time  thermoplastic  modulator  for  use  with  a coherent 
optical  processing  channel  was  accomplished.  Recording  bandwidths 
to  70  MHz;  spatial  bandwidth  of  30  c/mm;  dynamic  range  of  40  db; 
and,  high  optical  quality  recording  apertures  over  several  square 
inches  were  realized.  Optical  processing  for  the  spectrum  analy- 
zer application  with  either  rectangular  or  polar  input  data  re- 
cording formats  was  analyzed  and  the  unusual  properties  asso- 
ciated with  polar  data  were  clearly  established.  Fabrication 
and  coating  methods  for  a large  15"  diamater  rotatable  recording 
disc  and  5"  recording  plate  were  refined  to  give  thermoplastic 
recording  surface  of  high  thickness  uniformity. 

This  program  has  served  to  establish  the  feasibility  of  em- 
ploying a thermoplastic  modulator  as  the  means  for  high  quality 
and  high  density  data  insertion  into  a coherent  optical  processing 
channel  on  a real-time  basis,  thus  allowing  a more  practical  ap- 
proach to  utilization  of  the  excellent  properties  of  coherent  opti- 
cal processing. 

It  is  recommended  that  the  technology  base  established  in  this 
program  be  extended  in  the  areas  of  (1)  further  design  refinement 
for  increased  spatial  bandwidth  (=100  c/mm)  and  examination  of 
material  properties;  (2)  study  of  system  signal  response  for  in- 
creased dynamic  range  (=60  db) ; and  (3)  evaluation  of  additional 
real-time  optical  computer  application  areas  such  as  SAR  data  pro- 
cessing, pattern  recognition,  hybrid  optical- digital  systems  and 
special  purpose  display  devices.  The  thermoplastic  modulator 
should  also  be  considered  as  a device  usable  for  construction  of 
computer  generated  holographic  optical  elements  or  other  types 
of  thin  film  optical  components. 
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